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I \', ABSTRACT 

^^ i The intergalactic medium (IGM) is the dominant reservoir of baryons at all cosmic 

r) ' epochs. In this paper, we investigate the evolution of the IGM from z = 2 — >■ in 

• , (48/i~^Mpc)^, 110-million particle cosmological hydrodynamic simulations using three 

i-Q ■ prescriptions for galactic outflows. We focus on the evolution of IGM physical proper- 

»~^| ties, and how such properties are traced by Lya absorption as detectable using Hubble 's 

Q . Cosmic Origins Spectrograph (COS). Our results broadly confirm the canonical picture 

)—{ ' that most Lya absorbers arise from highly ionized gas tracing filamentary large-scale 

c/5 I structure. Growth of structure causes gas to move from the diffuse photoionized IGM 

Cu . into other cosmic phases, namely stars, cold and hot gas within galaxy halos, and 

the unbound and shock- heated warm- hot intergalactic medium (WHIM). By today, 

pq ! baryons are comparably divided between bound phases (35% in our favoured outflow 

^ ■ model), the diffuse IGM (41%), and the WHIM (24%). Here we (re)define the WHIM 

T-H ' as gas with overdensities lower than that in halos (p/p < 100 today) and temperatures 

OQ , T > 10^ K, to more closely align it with the "missing baryons" that are not easily de- 

■"nJ" ' tectable in emission or Lya absorption. Strong galactic outflows can have a noticeable 

CN I impact on the temperature of the IGM, though with our favoured momentum-driven 

ly-N . wind scalings they do not. When we (mildly) tune our assumed photoionizing back- 

f^ ' ground to match the observed evolution of the Lya mean flux decrement, we obtain line 

^^ I count evolution statistics that broadly agree with available (pre-COS) observations. 

^"^ ■ We predict a column density distribution slope of /(-/Vhi) oc iVj^j^™ for our favoured 

^ I wind model, in agreement with recent observational estimates, and it becomes shal- 

•rH . lower with redshift. Winds have a mostly minimal impact, but they do result in a 

^\f ' shallower column density slope and more strong lines. With improved statistics, the 

H I frequency of strong lines can be a valuable diagnostic of outflows, and the momentum- 

■ - - ■ driven wind model matches existing data significantly better than the two alternatives 

we consider. The relationship between column density and physical density broadens 
mildly from z = 2 — > 0, and evolves as p oc iV{j/'*10~°'^^^ for diffuse absorbers, con- 
sistent with previous studies. Linewidth distributions are quite sensitive to spectral 
resolution; COS should yield significantly broader lines than higher-resolution data. 
Thermal contributions to linewidths are typically subdominant, so linewidths only 
loosely reflect the temperature of the absorbing gas. This will hamper attempts to 
quantify the WHIM using broad Lya absorbers, though it may still be possible to do 
so statistically. Together, COS data and simulations such as these will provide key 
insights into the physical conditions of the dominant reservoir of baryons over the 
majority of cosmic time. 

Key words: galaxies: formation, large-scale structure, quasars: absorption lines, 
ultraviolet: general 
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1 INTRODUCTION 

Diffuse intergalactic hydrogen produces numerous weak 
redshffted Hi Lya (1216A) absorption features along 
lines of sight to distant bright objects suc h as quasars , 
the phenomenon k nown as the Lya forest (|Lvndsl Il97ll : 
ISargent et al.lll980l ). Advances in sensitive high-resolution 
spectroscopy and theoretical understanding of structure 
formation have led to the currently fa voured pa r adigm 
for the origin of the Lya fo re st (ICen et al.1 1994: 
Zhang. Anninos. fc NormanI 1 19951 : iMiralda-Escude et al.1 



19961 : iHernguist et al.lll996l ): it arises from highly pho toion- 
ized intergalactic hydrogen (jGunn fc PetersonI Il965l ) trac- 
ing uncoUapsed large-scale structure in the Cosmic Web. 
The general properties of the Lya forest are well-described 
by the Fluctuating Gunn-Peterson Approximation (FG PA ; 
[Weinberg. Katz. fc Hernauistlll997al : Icrofl et aLlll998l ). in 
which a tight relation between density and temperature in 
cosmologically-expanding photoionized gas (|Hui fc GnedinI 
1 19971 ) leads to a tight correlation between the Lya op- 
tical depth and the underlying mass density at each 
location along the line of sight. At z > 2, where 
the Lya transition falls into the observed-frame opti- 
cal, this approximation is quite goo d, and it can be 
used to infer the mea n baryon density (jRauch et al.l 119971 : 
[Weinberg et al.lll997bl 'l and the matter fluctuation spectrum. 



which in turn provides constraints on cosmological param- 



eters |Croft et al 



^ '1999, '200^; 'McDonald et al jboool . l2006l : 

IViel. Haehnelt. fc Springel 2004 : Viel et al.i20m 

In this pape r, a successor to iDave et al.l (|l999l ') and 
I Dave et al.l l|200ll '). we use cosmological hydrodynamic sim- 
ulations to model the IGM and the Lya forest at lower red- 
shifts. Here the Lya transition falls in the ultraviolet (UV), 
requiring more difficult space-based observations. Further- 
more, the physical description of the Lya forest departs 
more strongly from the FGPA, because filamentary struc- 
tures begin to grow so large as to heat a significant fraction 
of int ergalactic gas to well above photoionization tempera- 
tures IjDave et al. 1999i). This results in the emergence of the 
Warm-Hot Int ergalacti c Medium (WHIM; ICen fc Ostrikeil 
ll999l : lDaveeral .1119991 . 1200 ll : ICenfc Ostrikeil l2006l ). which 
contains a substantial fraction of all baryons today. Ad- 
ditionally, the evolving relationship between matter over- 
density and Lya column density is such that at low-z, a 
given overdensity produces a lower column density absorber: 
roughly speaking, a 10^*'^ cm~^ line at z ~ 2.5 is physicall y 
analogous to a 10^^ cm~^ line at 2 = (jDave et al.lll999h . 
This shift results in many fewer Lya forest absorbers to a 
given sensitivity limit, i.e. a much thinner forest of lines. 
All these factors make the low-z Lya forest more difficult 
to characterize both observationally and theoretically, an d 
not as useful for constraining cosmology (|Zhan et al.ll2005h . 
But the intergalactic medium (IGM) at low-z still contains 
the vast majority of baryons in the Universe, so understand- 
ing the \ow-z Lya forest and the WHIM remains a central 
challenge to assembling a complete picture of cosmic baryon 
evolution. 

One of the original Key Projects of Hubble Space Tele- 
scope was to understand the low-z Lya forest. The Faint 
Object Spectrograph (FOS), despite being unable to fully 
resolve Lya absorbers, provided critical insights into the 
nature of the \ow-z Lya forest. For example, even the earli- 



est observations showed that the number density evolution 
of Lya absor bers must slow dramatically between z > 2 
and z = (Bahcall et al.1 Il99ll : iMorris et al] Il99ll '). Ini- 
tially, it was suggested that \ow-z absorbers might repre- 
sent a different population than those at high-z, with the 
high-2 systems tracing the Cosmic Web and the low-z sys- 
tems arising in galaxy halos (jBahcall et al.lll996l ). This idea 
was supported by the apparently tight connecti on between 
colum n density and galaxy impa ct parameter IjChen et al.l 
Il998h . However. [Paveet al.l (|l999l ') showed that all these ob- 
servations can be fully accomodated within the structure 
formation model of the Lya forest. The change in num- 
ber density evolution can be explained as an ionization ef- 
fect (ITheuns. Leonard, fc Efstathioulll99a ). since the domi- 
nant sourc e of metagalactic flux (q uasars) diminishes rapidly 
at z < 2 JHaardt fc Madaul [Jooil ) . countering the effect of 
the declining mean density of baryons. The relationship be- 
tween Lya equivalent width and galaxy impact parameter 
arises because the IGM is denser around galaxies owing to 
matter clustering. Subsequent higher-resolution quasar spec- 
troscopy using the Space Telescope Imaging Spectrograph 
(STIS) has generally confirmed these interpretations, and 
STIS Lya fo rest data is well-des cribed by structure forma- 
tion models (jDave fc Tripdl200ll ). Hence the \aw-z Lya for- 
est is not fundamentally different than the high- 2: one, al- 
though the relationship of neutral hydrogen absorption to 
large-scale structure has shifted. 

The installation of the Cosmic Origins Spectrograph 
(COS) on Hubble should usher in a new age for the study of 
the low-z Lya forest. The dramatically increased sensitivity 
over STIS will enable a more robust characterization of the 
z < 1 IGM, and it will enable study of Lya absorbers at 
overdensities comparable to that detectable at high-z. COS 
will also enable tomographic probes of the IGM at trans- 
verse separations that will provide interesting c onstraints on 
the large-scale coherence of Lya absorbers (e.g. ICasev et al.l 
12008) . All these data will provide theoretical studies with 
new opportunities to understand the \ow-z IGM, and new 
challenges for models to meet. 

An advantage of low-z IGM work compared to high-z is 
that the galaxy population around \ow-z absorbers can be 
studied to much greater depth. The relation between the two 
provides unique insights into how the early and late stages 
of gravitational collapse are connected, and how the baryons 
have decoupled from the dark matter via shock heating and 
radiative cooling. The galaxy-absorber connection can be 
studied through the individual association of galaxies with 
absorbers or by statistical measures such as correlation func- 
tions. This is the subject of a subsequent paper in this series 
(KoUmeier et al., in preparation). Another paper in this se- 
ries (Oppenheimer et al., in preparation) will deal with metal 
line absorbers detectable using COS, and their relationship 
to galaxies and large-scale structure. 

In this paper, we employ state-of-the-art cosmological 
hydrodynamic simulations to study the z <1 IGM as traced 
by the Lya forest and WHIM. Where appropriate we include 
comparisons with existing data and predictions for COS. 
Our new simulations have 216 x more particl es and 80 x 
more v olu me compared to our pr evious studies in lDave et al.l 
(| 19991 ) and lDave fc Tripd l|200ll '). and they include more so- 
phisticated physical processe s compared to t hose w orks and 
more recent studies such as IPaschos et al.l l|2009f ) . For ex- 
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ample, our simulations now include metal-line cooling and 
a well-constrained heuristic model for galactic winds that 
enrich the IGM and regulate star formation in accord with 
observations. While w e will not discuss chemical en richment 
in this work (see e.g. lOppenheimer &: Davg|2008l ). the en- 
ergy input from such outflows can non-trivially impact the 
thermal conditions in t he IGM, particularly around galaxies. 
iTornatore et al.l ([2010') investigated the IGM in simulations 
similar to ours with the further addition of feedback from 
black holes, showing that it can significantly impact WHIM 
properties. 

In ij2]we describe our cosmological hydrodynamic sim- 
ulations and the generation and analysis of artificial spec- 
tra. In ^ we describe the physical conditions of the IGM 
from 2 = 2 — )■ 0, including how Lya absorbers trace large- 
scale structure and the emergence of WHIM gas. In !|4] we 
present statistics of Lya absorbers and their evolution. We 
also present comparisons to STIS data and predictions for 
COS, and study the impact of differences in spectral reso- 
lution and noise. In !j5]we discuss how Lya absorbers trace 
the physical conditions of the low-z IGM, including wide 
absorbers that trace WHIM gas. We summarize and discuss 
the implications of our results in ij6] 



2 SIMULATIONS 

2.1 The Code and Input Physics 

We employ our modified version (|Oppenheimer fc Davg 
I2OO8I') of the N-body+hydrodynamic code Gadget- 
2 ( Springe! l2005l l. Gadget-2 uses a tree-particle-mesh 
algorithm to compute gravitational forces on a set 
of particles, and an en tropy-conserving formulation 
ijSpringel fc H ernguist 2002) of Smoothed P article Hyd ro- 
dynamics (SPH; Gingold & Monaghan 1973; ILucvIIiQTTI ') to 
simulate pressure forces and shocks on the gas particles. 
We include radiative cooling from primordial gas as sum- 
ing ionization equilibrium (fo llowing Katz et al. 1996|) an d 
metal- line cooling based on Sutherland fc Dopital (1993). 
Star formation follows a 1 Schmidt (1959i ) Law calibrated 
to the Kennicutt (1998ab) relation; particles above a den- 
sity threshold where sub-particle Jeans fragmentation can 
occur are randomly selected to spawn a star with half 
their original gas mass. The interstellar medium (ISM) 
is modeled through an analytic subgrid recipe following 
iMcKee fc Ostriker |l977|). including en ergv returned from 
supernovae (ISpringel fc HernguistlbOOa ). 

Kinetic outflows are also included emanating from all 
galaxies that are forming stars. Gas particles eligible for 
star formation are randomly selected to be in an outflow. 
Outflowing particle have their velocity augmented by v^ in 
a direction given by vxa, where v and a are the comoving 
velocity and acceleration of the particle, respectively. The 
ratio of the probabilities to be in an outflow relative to that 
for forming into stars is given by r;, the mass loading fac- 
tor. Hydrodynamic forces on wind particles are "turned ofl" 
until the particle reaches one-tenth the threshold density 
for star formation, or a maximum time of 20 kpc/w^. This 
hydrodynamic delay allows wind particles to escape their 
host galactic disks, mimicking the effect of "chimneys" in 
the interstellar medium that are unresolvable in cosmolog- 
ical simulations but would allow the escape of gas in real 



galaxies. The choices for Vw and 77, including their possible 
dependence on galaxy or halo properties, define the "wind 
model" . 

We consider three wind models. The first has no 
winds, "nw". This model fails a number of observational 
tests, including o verproducing cosm ic star formation (e.g. 
iDave et al.|l200ll: iBalogh et al.ll200ll) and faihng to enrich 
the IGM (|Oppenheimer fc Davel l2006). but it serves as a 
baseline model to assess the impact of outflows on the dif- 
fuse IGM. 

The second wind model uses constant values oi rj = 
2 and Uw = 680 kms"'^ for all outflows from every 
galaxy; we call this our "constant wi nd" (cw) model. The 
cw mo d el is similar to t h at us ed in ISpringel fc Hernguistl 
(|2003bl ). ITornatore et all (|2010l ) (though their wind speed 
is somewhat lowe r), and the Ove rwhelmingly Large Sim- 
ulations (OWLS; IWiersma et al.l ^2009). This model uses 
all the supernova energy (assuming a Chabrier IMF) to 
drive kinetic winds. Doing so suppresses cosmic star for- 
mation to be broadly in agreement with observations, 
with th e added benefit of y i elding resolution-converged 
results (jSpringel fc Hernguistl l2003bi ) . In detail, it fails 
to enric h the IGM in agreement with high-z C iv ab- 
sorbers (|Oppenheimer fc Davell2006h or to enrich galaxies 
in ac cord with observations of the mass-m etallicity rela- 
tion (|Dave. Finlator. fc OppenheimeI1l2007^ ■ but it is still 
a physically reasonable scenario that captures some of the 
main effects implied by the observations. 

Our preferred model, which we refer to as 
"vzw", uses scalings expected for momentum-driven 
winds (jMurrav. Quatert. fc Thompson! 120051 ). though we 
note that such sca lings can be generated b y oth er phys- 
ical mechanisms (jDalla Vecchia fc Schavd 120081 ). Here, 
Vw = Scr-^/i — 1, and ri = cro/a, where a is the velocity 
dispersion of the host galaxy. We estimate a fr o m th e 
galaxy's stellar mass following IMo. Mao, fc Whitg ()l998h . 
using an on- the-fly friends-of-frien d s gal axy finder, as 
described in lOppenheimer fc Davg (|2008l ). The scaling 
factor /l is randomly chosen for each ejection event in 
the range 1.05 — 2, based on observ a tions of local star- 
bursts by iRupke. Veilleux fc Sandera (|2005r ). This choice 
also g ives outflow velocities consistent with that seen in 
z ~ 1 ()Weiner et aL 2009) and z ~ 2-3 (jSteidel et al.ll2004h 
star-forming galaxies. The ejection speeds are generally 
comparable to the galaxy escape velocities at low redshifts, 
but because wind particles are slowed by hydrodynamic 
interactions with surrounding gas, they typically do not 
escape their host galactic halos, at least in more mas- 
sive systems, ctq is a free parameter adjusted to broadly 
match the cosmic star formation history; here we choose 
ao — 150 kms^^. This yields 7; ~ 1 for L '^ L, galaxies at 
z ~ 2, which is in agreement with the constraints inferred 
from 2 ~ 2 galaxies by Erb (2008). 

The key features of our momentum-driven wind model 
are the higher mass-loading factors and lower ejection 
velocities in lower mass galaxies. This model has displayed 
repeated and often unique success at matching a range of 
cosmic star formation and enrichment data. This includes 
observations of IGM enrichment as seen i n 2 ~ 2 — 4 
Civ absorbers (|Opp enheimer fc Davell2006l. l2008l ). 2 ~ 
Ovi absorbers (|Oppcnlicimcr fc Dave 20 0^), and z ^ 6 
metal-line absorbers (|Oppenheimer. Dave, fc Finlatorl 
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l2009r ). Concurrently , it also matches the galaxy mass- 
metallicity relation (|Dave. Finlator. fc Oppenheimerl 120071 : 
iFinlator fc Davg 120081 ) , and the enrichment le vels seen in 
z = intragroup gas IjDave. Qppenheimer. fc Sivanandam 
l2008r ). It also suppresses star formation in 
agree ment with high-redshift lum inosity func- 
tions (|Dave. Finlator fc Oppenheimeill2006l ) and the cosmic 
evolu tion of UV luminosity at z ~ 4 — 7 (JBouwens et alj 
[20071). In contrast to the nw and cw models, the vzw model 
approximately reproduces the observed z = G stellar mass 
function of galaxies up to luminosities L ^ L^, though it still 
predicts excessive stellar masses for higher mass galaxies 
IjOppenheimer et al.l 120101 ). The numerical implementation 
of this model — with probabilistic, one-at-a-time particle 
ejection and the temporary turn-off of hydrodynamic forces 
— is only a heuristic representation of the underlying 
physics, and there are significant numerical and physical 
uncertainties associated with the hydrodynamic interactions 
between wind and ambient-halo particles and the mixing 
of outflows with the surrounding gas. Nonetheless, this 
heuristic model appears to distribute matter and energy 
into the IGM and regulate galaxy growth in broad accord 
with observations. 

While we leave a detailed discussion of metals to a forth- 
coming paper in this series, the simulations include a sophis- 
ticated chemical evolution model that follows four individ- 
ual species (C, O, Si, Fe) from three sources: Type II SNe, 
Type la SNe, and stellar mass loss from AGB stars. The for- 
mer instantaneously enriches star-forming particles, whose 
metallicity can then be carried into the IGM via outflows 
(descri bed below). Type la mod eling is based on the fit to 
data bv lScannapieco fc Bildstenl (|2005 ). with a prompt com- 
ponent tracing the star formation rate and a delayed com- 
ponent (with a 0.7 Gyr d elay) tracking stellar mass. Stellar 
mass loss is derived from iBruzual fc Gharlot (120031) popu- 
lation synthesis modeling assuming a lChabrierl ( 20031 ) IMF. 
Delayed feedback adds energy and metals to the three near- 
est gas particles; Type la's add lO^'^ ergs per SN, while AGB 
stars add no energy, only metals. We use lLimongi fc ChiefSl 
l|2005h yields for Type II SNe; yields for Type la and AGB 
stars come from various works (see lOppenheimer fc Davg 
|2008| . for details). 



2.2 Simulations 

The primary simula t ions used here are those described by 
lOppenheimer et al.l (|2010l ). using the three outflow mod- 
els outlined above. Each simulation contains 384^ gas and 
384"^ dark matter particles, in a random periodic volume of 
48/i~^Mpc (comoving) on a side with a gravitational soften- 
ing length (i.e. spatial resolution) of 2.5/i^^kpc (comoving, 
Plummer equivalent). We assume a cosmology concordant 
with WMAP-7 constraints (Ja,rgsik et al. 2010), speciflcally 
Qm = 0.28, ^A = 0.72, h = J^o/lOO kms~^ Mpc^^ = 0.7, 
Q.b = 0.046, n = 0.96, and erg = 0.82. This yields gas 
and dark matter particle masses of 3.56 x lO^Af© and 
1.81 X 10* M0, respectively. The initial conditions are gen- 
erated with an lEisenstein fc Hul (|l999l ) power spectrum at 
z = 129, in the linear regime, and evolved to 2 = 0; they 
are identical for all simulations. Our naming convention is 
r[boxsize]n[number of particles per side] [wind model], where 
the initial letter "r" indicates the particular choice of cos- 



mology above. Hence our models are r48n384nw, r48n384cw, 
and r48n384vzwl5, where the postscript 15 on vzw indicates 
our choice of ctq — 150 kms~^. 

To investigate the effects of the numerical resolution 
and box size, we also employ a simulation with 512'^ gas 
and 512^ dark matter particles, in a random periodic vol- 
ume of 96/i~^Mpc (comoving) on a side with a gravitational 
softening length of 3.75/i~^kpc (r96n512vzwl5). This simu- 
lation is run with the same momentum-driven wind model 
as our 48/i^^kpc case, and we will refer to it as "vzw-96". 
Compared to vzw-48, it has an eight times larger simula- 
tion volume but 3.375 times worse mass resolution. Hence it 
tests both sensitivity to box size and numerical resolution, 
although the dynamic range in either sense is not dramati- 
cally large. 



2.3 Generating Spectra 

Spectra are generated using our spectral generation code 
specexbin, which casts lines of sight at an angle to the sim- 
ulation axes, allowing continuous and non-repeating lines of 
sight to be generated through the periodic simulation vol- 
umes. We generate 70 such continuous spectra from z — 2 
to z — 0. specexbin calculates optical depths in Hi and a 
variety of metal ions; here we will only be concerned with 
Lya. For line statistics we will typically consider three red- 
shift ranges, namely z = — 0.2 (which we will call "z = 0"), 
2 = 0.9 - 1.1 ("2 = 1") and 2 = 1.8 - 2 ("2 = 2"). In each 
interval, the total path length for the 70 spectra is Az = 14. 
This is comparable to the path length to be obtained by, 
e.g., the program of the COS GTO team. 

The procedure for calcu lating optical depths is de- 
scribed by Oppcnh eimer fc D ave (2008). In brief, the algo- 
rithm is to (i) calculate the neutral fraction of each gas 
particle based on the density, temperature, and the as- 
sumed photoionizing flux (discussed below); (ii) smooth 
the neutral hydrogen density onto randomly-chosen lines 
of sight in real space; (iii) calculate Hi optical depths as- 
suming ionizatio n equilibrium and optically-thin gas us- 
ing CLOUDY96 (|Ferland et al.lll998l ): and (iv) use the Hi- 
weighted peculiar velocity as well as Hubble flow to rebin 
the real space optical depths into redshift space. This pro- 
duces a table of optical depths versus redshift along each 
line of sight. Concurrently, for each ionic absorption species 
specexbin outputs the optical depth-weighted temperature, 
density, and metallicity, which are used to quantify the phys- 
ical conditions in each identified absorber. For lines of sight 
that parallel box axes, spece xbin produces the sa me Lya 
absorptio n spectra as tips-^ jHernguist et al.lll996l ). As re- 
ported bv lPeeples et al.l (|2010h . the Lya flux power spectra 
computed from tipsy spectra at 2 = 3 are indistinguishable 
from th ose produced by t he independent spectral extraction 
code of iLidz et al.l l|2009l ). confirming the robustness of the 
procedure to the details of numerical implementation. 

Extracted spectra are converted from optical depths to 
fluxes, resampled to a pixel scale of 6 kms~^, and convolved 
with the COS line spread functiorjj (LSF). Speciflcally, we 
use the FUV G130M LSF at 1450A. The LSF varies only 



^ http : //www-hpcc ■ astro ■ Washington ■ edu/tools/t ipsy/tipsy ■ html | 



^ http://www.stsci.edu/hst/cos/performance/spectraLresoIutionl 
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Figure 1. Top: Mean flux decrement (D) evolutio n from z = 2 —^ 
in o ur three simulations. Data points are from lKirkman et al.l 
1120071 ) . The numbers in parantheses are the values of /^ by which 
the optical depths are multiplied in each simulation to obtain 
the agreement with data as shown. Bottom: Fhi evolution that 
matches the mean flux for the vzw (green) and nw (blue) runs. 
The cw run (not shown) is identical to nw. These are simply the 
iHaardt fc Madaul 1 120011 ) values divi ded by fr- Data p oints show 
proximity effect measurements from lScott et al.l 1 120021) . 



very slightly with wavelength in the FUV channel, but is 
substantially different (and narrower) in the NUV channel. 
Nevertheless, we keep the LSF fixed at all wavelengths so as 
to avoid introducing any artificial evolution from the varying 
LSF. Generally, the FUV LSF is somewhat broader than 
had been originally expected, with its central peak having 
an equivalent la Gaussian width of ~ 17 kms"^, and it 
also has substantial non-Gaussian wings. Finally, we add 
Gaussian random noise to the simulated spectra assuming a 
signal to noise ratio of S/N=30 per pixel. This is comparable 
to the highest quality data that will be obtained with COS. 
The shape and amplitude of the metagalactic pho- 
toionizing flux are still not well constrained, particularly at 
lower redshifts. We assume that the flux is spatially-uniform, 
which is probably a good approximation in the diffuse and 
optically thin Ly g forest. We take the spe ctral shape from 
that predicted by iHaardt fc Madaul l|200ll ) using the popu- 
lation of observed quasars and star-forming galaxies (their 



"QlGOl" model). This model assumes a 10% escape fraction 
of ionizing radiation from gala xies and is broadly consistent 
with available constraints (e.g. IScott et al.ll2002l ). 

We adjust the amplitude of the metagalactic flux to 
match the observed evolution of the mean flux decrement 
(D) in the Lya forest. A modest post-processing adjust- 
ment recovers essentially the identical answer to running 
the entire simulation with a different ionizing background, 
since the dynamics of Lyg forest gas are n egligibly affected 
by photoionization (jWeinberg et al.lll998l ). The amount of 
adjustment varies slightly among our simulations, since the 
outflows have a small but non-negligible impact on the Lya 
forest. In practice, we multiply all the optical depths by a 
factor /r before we convert to flux, then measure D, and 
iteratively adjust fr until we obtain a good match to D{z). 
In the optically-thin approximation (excellent in this case), 
this procedure is equivalent to multiplying the strength of 
the ionizing background by 1//t. While in principle we could 
have f-r be a function of redshift, this was not seen to be nec- 
essary to within the accuracy of available data, and so we 
simply determine a constant /t- individually for each simu- 
lation. 

The resulting evolution of the mean flux decrement 
in the Lya forest is shown in the top panel of Figure [T] 
The values of /t are indicated in the parentheses within 
the lege nd; they are close to u nity, indicating that the 
original IHaardt fc Madaru (|200ll ) metagalactic flux ampli- 
tude is viable to within current observational uncertainties. 
In all cases, th e resulting D(z ) is a reasonable match to 
observations by iKirkman et al.l (|2007l ). who measured the 
flux decrement evolution using data compiled from Hubble^s 
Faint Object Spectrograph (FOS) at z < 1.5 and the Keck 
telescope's High Resolution Echelle Spectrograph (HIRES) 
a,t z > 1.5. The simulations appear to have a somewhat 
higher flux decrement at 1 < z < 1.5, but in this regime the 
low resolution of FOS and the increased blending owing to a 
thicker Lya forest (as compared to 2 ~ 0) are more likely to 
result in continuum placement errors that would cause the 
flux decrement to be biased low. Hence we do not consider 
these discrepancies serious as yet, though they are worth 
revisiting when better data becomes available from COS. 

The bottom panel of Figure [1] shows the evolution 
of the Hi photoionization rate Fhi that yields the mean 
flux evolution in the top panel. Since we multiply the 
IHaardt fc Madaul (J200ll ) model by a constant factor, the 
shapes of the curves are identical to their model. The am- 
plitude is multiplied by fr^, so vzw is slightly higher than 
nw; cw is identical to nw since it has the same fr- All the 
models are consistent with the availabl e observational con- 
straints as they are fairly uncertain (see lDave fc TripDll200ll . 
for more constrain ts); here we sho w proximity effect mea- 
surements from S cott et al.l (|2002l ). If we take the z — lA 
measurement of D at face value and modify Fhi (2) to match 
it, this would result in an increase in Fhi by 0.2 dex at 
z — 1.4, and a subsequent recovery back to the original 
track by 2: = 1.6. This would be an odd detour in its evo- 
lution, implying a peak in photoionizing flux at 2 ~ 1.5 
alo ng with a seconda r y pea k at 2 ~ 2 — 2.5 as in the origi- 
nal [Haaxdt_&_^^adau| (200l|) evolution. In this case the peak 
epoch of Ji, would no longer correspond to the peak epoch of 
high-luminosity quasars, although lower luminosity quasars 
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Figure 2. Example of a section of noise-free simulated spectra (A2 = 1), with an expanded region shown on the right, comparing three of 
our outflow models. Top panel shows Lya absorption, middle shows density relative to the cosmic mean, and bottom shows temperature. 
The color code is nw (blue) plotted on bottom, cw (red) above, and vzw (green) on top; e.g. regions where only green shows means the 
others are identical to within the display resolution. The density structure is similar in all models, with the vzw having slightly more 
absorption overall. The cw model clearly produces a warmer IGM and the most deviations in optical depth; the other two show smaller 
differences. 



peaking at lower redshifts could provide a siRnifican t contri- 
bution (jShankar. Weinberg. Miralda-Escude II2009I ). 

It is interesting that the constant wind and no wind 
cases have identical values of /t, while the vzw case is lower. 
In the simplest view, winds tend to add energy to the IGM, 
resulting in less Lya absorption, and hence one needs smaller 
optical depths and a lower /r. In detail however, because D 
is dominated by marginally saturated absorbers, the relevant 
issue is whether the winds are heating gas giving rise to those 
absorbers in particular. It turns out that vzw, owing to its 
smaller wind speeds from typical galaxies, tends to deposit 
energy closer to galaxies where strong absorbers arise, while 
cw tends to deposit energy far away in more diffuse gas, as 
we will discuss in 33] Hence we believe that although cw 
adds more energy to the IGM overall, this does not impact 
D nearly as much as in vzw. 

Figure [5] shows an example spectrum from z — 0.3 — 0.4 
along an identical line of sight (LOS) through our three wind 
simulations. The Lya forest is fairly sparse at these epochs, 
but many very weak lines exist that could be detected given 
sufficient observational capabilities. The mean absorption is 
similar by construction, but there are distinct variations ev- 
ident particularly between the cw case and the other two. 
The expanded segment in the right panels highlights a region 
that shows significant differences among all three wind mod- 
els. The density structure is mostly independent of outflows, 
though small differences are evident, but the IGM temper- 
ature can be significantly impacted by wind energy. The 
constant wind case shows a higher temperature almost ev- 
erywhere along the LOS, while the temperatures for vzw and 
nw are fairly similar. However, the fact that optical depths 
are not very sensitive to temperature makes these differences 
difficult to pick out in the observed spectra. Overall, winds 
do not have a strong impact on the physics or absorption 
properties along typical lines of sight, as we will demonstrate 
more quantitatively with Lya absorber statistics. 



3 IGM PHYSICAL CONDITIONS 

3.1 Cosmic Phase Diagram 

The left panels of Figure [3] shows cosmic phase diagrams 
of baryons, i.e. overdensity versus temperature, in our fidu- 
cial vzw run at redshifts from 2—^0. The trends shown 
are familiar. At the lowest densities there is a tight density- 
temperature relation set by a competition between photoion- 
ization heating and adiabatic cooling owing to Hubble ex- 
pansion. At roughly the cosmic mean density, shock heating 
begins on collapsing filamentary structures, creating a plume 
of hotter gas. At higher densities, radiative cooling becomes 
effective, creating the spur of dense gas that fuels star for- 
mation in galaxies. The other wind models look qualitatively 
similar. 

Two perhaps unfamiliar features are the shelf of gas 
sitting at T « 10''K, which is metal-enriched gas outside 
of star formation regions (we truncate all cooling at lO^K), 
and the upwards-curving plume of the densest gas, which is 
a consequen ce of the two-phase interstel lar medium imple- 
mentation of lSpringel fc Hernguistl ([2003). The temperature 
in the two-phase regime is given by a mass-weighted aver- 
age of the subgrid hot phase, assumed to be at 10* K, and 
the mass-dominant cold phase at 1000 K. Neither of these 
features are physical; they simply represent computational 
conveniences. 

From redshift 2 to 0, we see the growth of galaxies and 
large-scale structure playing out in the cosmic phase dia- 
gram. The denser phases become more heavily populated, 
as does the plume of shock-heated WHIM gas, though the 
hotter phases are still clearly populated even at z = 2. The 
hottest gas also grows most substantially owing to the for- 
mation of large potential wells at late times. The diffuse gas 
extends into lower overdensity regions at late times, as voids 
become more empty, and its temperature drops with time 
as the physical density and photoionization rates go down. 
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Figure 3. Left panels: Cosmic phase space diagram for gas particles in tlie vzw run at z = 2, 1, (2000 randomly-subsampled gas 
particles). We divide the baryons into four categories: Diffuse, Condensed, Hot Halo, and WHIM. This is done using a temperature 
cut at 10^ K to separate the former two from the latter two, and a density cut given by equation \T\ to separate the bound baryons 
(condensed and hot halo gas) from the intergalactic baryons (diffuse and WHIM). Temperature histograms are along the j/-axis for each 
phase and the total (black), with median values indicated by the tick marks. The percentage of mass in each phase is listed in each panel; 
additionally, the fraction in stars is shown in magenta. Right panels: Same, for the other simulations (nw, cw, and vzw-96), all at ^ = 0. 



The overall cosmic gas temperature distribution becomes 
broader with time, and shifts to slightly higher tempera- 
tures, driven by the increasing fraction in WHIM and hot 
halo gas. These tren ds have all been n oted and described in 
previous works (e.g. iDave et al.iri999l ). 



3.2 Redefining Cosmic Baryon Phases 

In Figure [3] we have also subdivided our baryons into four 
cosmic phases. The su bdivisions are slig hly different than 
conventional ones (e.g. iDave et al.l 120011 ) . in which diffuse 
baryons are at T < 10^ K, WHIM gas is 10^ < T < lO'^ K, 
hot gas is T > 10^ K, and condensed baryons are at a fairly 
high overdensity. Our new divisions are intended to be more 
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physically-motivated, particularly in the case of the WHIM, 
though in practice they do not represent dramatic changes. 
The main difference is that we divide the phases accord- 
ing to an overdensity threshold. This threshold is intended 
to reflect the overdensity (5 = p/p — 1) at the boundary of 
a virialized halo. To determine this, we note that fo r rea- 
sonable choices of the iNavarro. Frenk. fc White! l|l996r ) con- 
centration parameter c, the density at the virial radius is 
roug hly one-third the mean e nclosed densitjij. We thus ob- 
tain (|Kitavama fc SutdlToQl ) 



where 



fn 



5th = 67r^(H- 0.4093(l//n - 1) 



fi„(l + z)3 + (1 _ n™ - nA)(l + zY + Qa ' 



(1) 



(2) 



5th serves as our division between "bound" phases (gas in 
galaxies and halos) and "intergalactic" phases (WHIM and 
diffuse gas). The value of 1 -I- Sth evolves from Bvr^ « 60 at 
1 -I- « > (nA/fim)^^^ to « 120 at z = 0. This boundary is 
not completely sharp, as some gas that has higher overden- 
sities is not yet bound, and owing to the elliptical nature of 
halos some bound gas can have lower overdensities than Sth- 
We checked thi s definition by run ning a spherical overden- 
sity halo finder l|Kere8 et al.ll2005i l , and checking that equa- 
tion U]) faithfully and robustly separated gas within halos 
from gas outside. It indeed did so, as unbound gas was never 
more than a couple tenths of a dex in density into the bound 
region, and vice versa. 

We further separate the "hot" phases (WHIM and hot 
halo gas) from "cool phases" (diffuse and condensed gas) 
by a temperature threshold of Tth = 10^ K, the same value 
used in earlier definitions. At high densities, all gas that is 
star-forming is included in the condensed phase, although 
the subgrid two-phase model creates a mean gas particle 
temperature that can exceed T > Tth. At IGM densities, 
this represents the temperature above which the H I neutral 
fraction starts to drop dramatically, making Lya absorption 
a poor tracer of its baryonic content (though as we will see 
later, wide Lya absorbers can still trace this gas). Further- 
more, since photoionization cannot easily heat gas to these 
temperatures, it demarcates gas that has been shock-heated 
by large-scale structure. From Figure [S] one can see that 
some gas at T = 10^ — 10^ K is also shock-heated above 
the photoionization locus, but the locus itself extends to 
T Ri lO'^ K, at least at z ^ 1. 

With these divisions, we have four phases as follows: 

• Diffuse (r< Tth, 5 < (5th), 

• WHIM (T> Tth, 5 < 5th), 

• Hot halo (T > Tth, 5 > 5th), 

• Condensed (T < Tth, 5 > 5th), 

with the qualification that two-phase star-forming gas is 
"condensed" regardless of temperature. The most significant 
departure from previous definitions is that of the WHIM gas. 
Unlike the common definition in which 10^ — lO'^ K gas is 



WHIM regardless of density, the WHIM gas in our defini- 
tion is truly intergalactic, in the sense that it is not within 
the boundaries of a galaxy or its halo. Also, since hot gas is 
most commonly detected in X-ray emission, which roughly 
scales as oc p'^T^'^ (neglecting metal line emission for the 
sake of this argument), this lower-density gas is more likely 
to be "missing" from current observational censuses. For in- 
stance, our definition removes hot gas in galaxy groups at 
T < 1 keV from being part of the WHI M; such gas can be de- 
tected in soft X-rays (|Mulchaevll200Gl ) because its high den- 
sity (relative to truly intergalactic gas) leads to significant 
emission, albeit at faint levels. Hence, our revised definition 
more closely reflects the idea that the WHIM is the reposi- 
tory of the "missing" cosmic baryons. Eliminating the upper 
limit of T = 10^ K for the WHIM has minimal impact, since 
very little gas with 5 < 5th has such high temperatures. 



3 In detail, p(r2oo)/p 



3 (l + c)^(in(l + c)-c/(l + c)) 



wliicli is ap- 



proximately 1/3 for low values of concentrations corresponding 
to objects just undergoing collapse. 



3.3 Winds and the Phase Diagram 

The right panels of Figure |3]show the 2 = cosmic phase di- 
agram for our other simulations: nw, cw, and vzw-96. Com- 
paring vzw-96 to vzw, we see some subtle differences partic- 
ularly in the high-temperature gas. The hot halo gas extends 
to slightly higher temperatures, owing to larger halos being 
able to form in the 96/i~'^Mpc volume; this causes a slightly 
higher median temperature for this phase, from 10^'*^ K 
to lO'^'^^ K. Similarly, the WHIM gas extends to slightly 
higher temperatures, causing a comparable increase in me- 
dian WHIM temperature. The cooler phases show no ap- 
preciable physical differences. There is also a feature in the 
condensed phase at p/p « 10® corresponding to a slightly 
different implementation of the two-phase cutoff in this run, 
but this does not affect any of our results. 

The no-wind case also shows some differences with re- 
spect to vzw. Although the distribution of WHIM gas in 
phase space is similar, the median temperature increases 
slightly from 10^'^'' K to lO'^'*^ K, just as in the vzw-96 
case but for a different reason. Here, the main effect is 
that w ith no outflows, the IGM is almost completely unen- 
riched (|Oppenheimer fc DavgbOOST ). Hence there is no metal 
cooling in this regime, which yields somewhat higher tem- 
peratures. For the same reason, the condensed phase shows 
almost no particles at the 10'' K metal cooling floor. 

The constant wind model noticeably pushes WHIM gas 
into a lower-density regime compared to vzw, but its median 
temperat ure is lower (10^'®^ K) since t he diffuse IGM is more 
enriched (|Oppenheimer fc Davell2006h . In contrast, the hot 
phase temperature is higher by « 0.15 dex, reflecting the sig- 
nificant IGM heating that occurs when winds from all galax- 
ies, even small ones, are emanating at 680 k ms~'. These re- 
s ults a re qualitatively similar to those from iTornatore et al.l 
C201CI ) who explore a similar wind model, and they further 
note that black hole feedback (as they implement it) causes 
even more heating of particularly low-density gas. It is also 
similar to results from mesh hydrodynamic simulations by 
ICen fc Ostrikeri (|2006l ). who found that winds increase the 
mass fraction in the WHIM by 20%. Hence, it is possible 
to non-trivially impact the cosmic distribution of baryons in 
phase space with plausible outflows. 
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Figure 4. Top left: Evolution from z = 4 — > of baryon fraction 
in various cosmic phases, including stars, in our fiducial vzw run. 
Next 5 panels: Evolution of baryon fraction in stars (upper right), 
and then the various baryonic phases in order of their location 
in p — T space: WHIM (middle left), hot halo (middle right), 
diffuse (lower left), and condensed (lower right), in our three wind 
models (green, blue, red lines), along with the larger- volume vzw- 
96 simulation (cyan). 



3.4 Baryon Phase Evolution 

We now quantify the trends in the cosmic phase diagram by 
examining the evolution of baryons in each phase amongst 
our various simulations. As is apparent from Figure [3l 
baryons progressively move from the diffuse phase into 
WHIM and hot halo gas, owing to the hierarchical growth 
of large potential wells. The condensed phase gas is supplied 
predominantly fro m cold mode accre tion but also some cool- 
ing from halo gas (jKeres et aLlbOOg ). and this provides the 
fuel to grow the cosmic stellar content. By z = 0, baryons 
are comparably divided between diffus e, WHIM, a nd th e 
various bound phases, as was found in iDave et al.l (|200lh . 
However, the new definition of the WHIM yields somewhat 
less WHIM gas and more diffuse gas, tilting the balance to- 
wards diffuse gas overall. The precise balance between the 
diffuse and WHIM phases is also sensitive to the adopted 
value of Tth , as one can see from Figure (3] 

Figure |4l top left panel, shows the fraction of baryons 
in each cosmic phase from 2 = 4 — > in our vzw run. We 
additionally show the evolution of baryons in stars. Quali- 
tatively, the behavior is very similar to that seen in all cos- 



mological hydrodynamic simulations (jDave et al.ll200ll ). as 
discussed above. The next five panels of Figure |4] show a 
comparison of individual phase evolution for our three wind 
prescriptions, plus our large-volume vzw-96 run. The upper 
right panel shows the evolution of the global stellar con- 
tent, i.e. the integra l of cosm ic average star formation rate 
density vs. redshift (|Madau| [l998;) , in each wind run. Note 
that these values account for stellar evolution which returns 
material back into the gas phase. The no-wind case pro- 
duces by far the most stars, which is the well-known result 
that with out strong feedba ck, there is a severe overcooling 
problem IjDave et al.l I2OOII ). A stellar baryon fraction to- 
day of ~ 2 0% far exceeds current observational values of 
< 10% (e.g. lBell et"ai1l2003l ). Though we don't show it here, 
thi s case is also poorly resolution-converged, as pointed out 
by iBalogh et alj (l200l|), and hence the problem would be 
worse at higher resolution. 

By contrast, both wind models do a good job of sup- 
pressing star formation. Their globally averaged star for- 
mation rates are almost identical down to z ~ 2, which 
is by design as the wind prescriptions are tuned to repro- 
duce the observed star formation rates up to this epoch. 
At lower redshifts, there start to be greater differences. Be- 
cause the vzw model has more recycling of wind material 
back into galaxies, it provides more fuel for late-time star 
formation ( Oppenheimer et al.ll201Gl ). None of our simula- 
tions includes a prescription (e.g., AGN feedback) to trun- 
cate star formation in the most massive systems as ob- 
served, which would lower the predicted stellar fractions 
particularly at late times. The resolution dependence of 
global star formati on is mitigated by the prese nce of out- 
flows, as shown bv iSpringel fc Hernguistl (|2003bj) for the cw 
case, and as shown here by comparing vzw and vzw-96. 
Both wind models are generally consistent with the cur- 
rently favo ured values of the stellar baryon co ntent today 
of 4 - 8% (iBaldry. Glazebrook. fc Driver Il2008h . 

The middle left panel shows WHIM evolution, which 
displays significant variations, particularly between models 
with and without galactic winds. The modest wind speeds of 
the vzw model do not heat much gas to WHIM temperatures 
compared to the purely gravitational (nw) case, but there 
is still some heating evident at high-«, which diminishes to 
z — 0. The cw case shows a significant amount of WHIM gas 
even at z = 4, and the WHIM fraction increases faster than 
in the other cases down to z ~ 1. At this point the supply 
of wind energy abates, and the W HIM fraction levels off 
and actually falls slightly to z = 0. iTornatore et al.l ([20101) 
found less dramatic effect on the WHIM fraction with similar 
simulations, although the energy input from their winds is 
half of ours. These differences highlight an important point: 
galactic winds can in principle have a strong effect on the 
evolution of WHIM gas. The winds in our favoured model 
do not substantially affect the WHIM, but faster winds in 
plausible variant models can do so. Measuring the evolution 
of the global WHIM content could provide constraints on the 
large-scale thermal impact of outfiows. Such measurements 
will be challenging given the difficulty of quantifying the 
WHIM even at the present epoch, but the redshifting of 
some UV tracer lines to longer wavelengths could potentially 
facilitate such measurements. 

The middle right panel shows hot bound gas evolution. 
In general this is dominated by gravitational heating from 
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large potential wells, as indicated by the no-wind model. The 
vzw case produces more hot bound gas at late times, as it 
removes material from galaxies and places it into their sur- 
rounding hot halos. Since the wind speed scales with escape 
velocity, large galaxies can have outflow velocities that ex- 
ceed the wind speed of the cw model, which provides greater 
heating in their large halos. The cw model removes a signifi- 
cant amount of hot gas at early times because its high-speed 
winds can escape smaller potential wells, but at late times 
the potential wells become deep enough to retain the gas, 
and its hot phase rises more quickly. Still, even by 2: = 
this model does not produce as much hot bound gas as the 
vzw model. The vzw-96 and vzw simulations are very sim- 
ilar again, showing that simulation volume effects are not 
important even for large potential wells in terms of their 
global mass content. 

In the lower left panel, the diffuse phase shows simi- 
lar trends in all the models. The vzw winds have little im- 
pact on the diffuse gas content, as their typically modest 
velocities do not produce substantial IGM heating. The cw 
winds lower the diffuse gas fraction at high-z by shock heat- 
ing some IGM gas above 10^ K. This heating produces no- 
ticeable broadening of high-z metal lines, leading to worse 
agreement (relative to vzw) with observed C iv line widths 
l|Oppenheimer fc Davg 120061 ). At late times, heating of the 
low density IGM declines because of the reduced global star 
formation and the trapping of the constant-u„ winds by 
deeper potential wells, causing the diffuse phase fraction to 
return to the no-wind value by z ~ 0.5. The box size has 
a negligibly small impact on this phase — the cyan vzw-96 
curve is hidden by the overlaid green vzw curve at all but 
the highest redshifts. 

The bottom right panel shows the condensed phase evo- 
lution. The overall trend for this phase is to increase at 
early epochs as material accumulates in halos, then decrease 
towards later times as galaxies convert gas to stars more 
rapidly than they accrete fresh gas. The fraction of con- 
densed gas is < 10% at all times and < 6% today, in all 
the models. The cw model shows a fixed offset to lower con- 
densed gas fractions relative to the no-wind model, since it 
expels a constant amount of material out of the star-forming 
regions, and most of this material ends up outside of halos 
owing to the large wind speed. The vzw model shows quali- 
tatively different behavior from the others, with less redshift 
evolution. Mostly this difference refiects the late-time reac- 
cretion of wind material, more p rominent in the vzw mode l 
because of lower wind speeds IjOppenheimer et al.l |2010^ . 
We note the qualitative similarity of the evolution of this 
phase (which contains most of the cosmic neutral gas) with 
the obser ved evolution of co smic H I gas from Damped Lya 
systems (| Wolfe et al.l |2005| ). being roughly constant from 
~ 3 — >■ 0.5 and then dropping to z = Q. 

Returning to Figure O the histograms along the verti- 
cal axis show the temperature distribution of baryons at 
z = 2,1,0. The overall trend follows the movement of 
baryons within the cosmic phase diagram, namely that the 
median temperature increases with time. In detail, the dif- 
fuse phase actually moves towards lower temperatures be- 
cause of the changing balance between photionization heat- 
ing and adiabatic cooling. The WHIM gas extends up to 
higher temperatures at low-z, and the hot bound gas also 



becomes hotter on average as larger potential wells form. 
The median cosmic gas temperature today is « 10^ K. 

Overall, the cosmic evolution of baryon phases refiects 
the hydrodynamic and radiative processes associated with 
the growth of structure. Winds have a significant impact, 
particularly on the stellar phase (by design) and on the 
WHIM phase (if the winds are strong enough). They also 
impact the bound phases (hot halo and condensed), since 
these lie close to galaxies where winds play a more impor- 
tant dynamical role. The diffuse phase is mostly unaffected 
by winds. The complex dynamical interplay between out- 
flows and hierarchical growth is an important factor in the 
cosmic history of baryons, and improved observational con- 
straints and theoretical predictions will be required to fully 
characterise it. 



4 Lya ABSORBER OBSERVABLES 

We now investigate the statistics of Lya absorbers in our 
simulated COS-like spectra, and present comparisons when 
appropriate to existing pre-COS data. We consider the 
canonical statistics of the column density and linewidth dis- 
tributions, as well as the evolution of the number density 
of lines. To obtain line statistics, we fit o ur simulated spec - 
tra with Voigt profiles using AuTOVP l|Dave et al.lll997l ). 
yielding an Hi column density (A^hi), Doppler line width 
(b), and redshift for each absorber. Owing to the sparse- 
ness of the low-z Lya forest, the fits are generally unam- 
biguous, and tests have shown that AuTOVP obtains line 
parameters very similar to that from other Voigt profile fit- 
ters (e.g. VPFIT) when the lines are unsaturated. We use 
the continuum directly provided by the simulations, as we 
assume that observed spectra can be fit with an accurate 
continuum in the relatively sparse low- 2; Lya forest. We also 
examine the impact of varying the assumed signal-to-noise 
ratio and spectral resolution of our simulated spectra. Note 
that AutoVP provides "raw" &-parameters uncorrected for 
instrumental resolution. In our COS-resolution spectra, with 
a Gaussian width of 17 kms~^, we therefore expect essen- 
tially no lines with b < 17^2 ^ 24 kms~^. Since it turns 
out that the vzw-96 absorber statistics are indistinguishable 
from the vzw case (as anticipated by the close agreement in 
IGM physical properties) , we do not show this model herein. 



4.1 Column Density Distributions 

Figure [5] shows column density (left panels) and 6-parameter 
(right panels) distributions for absorbers at z — 0,1,2. As 
mentioned previously, these three bins correspond to all lines 
within redshift ranges z = — 0.2, z — 0.9 — 1.1, and 
z = 1.8 — 2, respectively, identified with AuTOVP in our 
70 random lines of sight. Results are shown for all three of 
our outflow models. The column density distribution (CDD) 
/(A^Hi) = cPN/dNuidz has been multiplied by A^hi to re- 
move some of the steep power-law scaling (discussed below) , 
which enhances the visibility of the differences among mod- 
els. It also results in a dimensionless number that reflects 
the relative number of lines per unit redshift at a given A^hi. 
The CDD shows the usual power-law behaviour, above 
the completeness limit at TVhi ~ lO^'^cm"'^ that results 
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Figure 5. Column density (left) and b-parameter (riglit) distributions at z = (top), z = 1 (middle), and z = 2 (bottom) for our 
three simulations. The CDD have been multiplied by N m improve visibility of the differences between models. The thick cyan line 
shows a fit to data from iPenton. Stocke. &: Shulll 1 12004) at z ss 0, namely log/(A'^Hl) = 10.3 — 1.65 log A^jjli over their fitted range of 
12.3 ^ logTVui ^ 14.5; dotted lines on either side indicate thei r la q uote d uncertainty of ± 1.0 in log/(AfHi). Thick black and magenta 
lines show best fits for CDDs from STfS data from lLehner et al.l 1 120071 ) and lDave et al.l l l200M . respectively. The same lines are reproduced 
at the other redshifts for reference, showing that there is modest evolution in the CDD amplitude from z = 2 — >■ 0. The turnover at 
A^HI < lO^^cm"^ arises from incompleteness in recovering lines at the assumed S/N and resolution; see t|4.4l Linewidth distributions are 



shown for all lines with TVjji > 10 cm 
curve indicate median b values. 



our complete sample at COS resolution and an assumed S/N=30. Tick marks above each 



from our adopted S/N and resolution (see ^4.4p . Charac- 
terizing the column density distribution as a power law 



with f{Nj 



m) oc 



NZ 



and fitting between 10" < A^^hi < 



10 "'cm , we obtain values oi P — — 1. 



-1.70, -1.80 for 



the nw, vzw, and cw models at z = 0, respectively, with 
an uncertainty of approximately ±0.1 on each. Hence, all 
the slopes are formally consistent with each other, although 



given that the exact same LOS are being analyzed, one 
can conclude that the nw model produces a steeper slope 
than the wind models, and that the vzw model produces 
the shallowest slope. We note that all the slopes are shal- 
lo wer than the slope — 2.04±0.23 obtained in the simulation 
bv lDave et al.l (|200in (though within formal uncertainties), 
which was most analogous to the no-wind case here. We will 
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show in §43]that higher spectral resolution tends to produce 
a somewh at steeper slope, so t his could explain part of the 
difference JPasch os et all 12009*) found a slope of -1.84±0.02 
using a fixed mesh simulation with 75 kpc resolution and no 
galactic outflows, which basically agrees with our no-wind 
case. 

These are several pre-COS set of low-z Lyg forest data 
for co mparison. One sample, from IPenton. Stocke. &: Shulj 
l|2004 ). contains 187 somewhat heterogeneous Lyo ab- 
sorbers, mostly at 2 < 0.1, from Hubble/STIS and Goddard 
High Resolution Spectrograph data. A more uniform and 
higher-resol ution sample o f 341 absorbers out to z < 0.4 was 
obtained bv lLehner et al.l (|2007r ) using only STIS data, and 
included Far Ultraviolet Space Explorer data to help con- 
strain the properties of s aturated Lya abs orbers. We also 
show earlier results from iDave et al.l l|200ll ) , who analyzed 
two STIS spectra using AuTOVP. The fits to all their re- 
sulting CDDs are shown a s the thick lines in the upper left 
panel of Figure[5l black for lLehner et al.l ((20071), magenta for 
I Dave et al.l l|200ll 'l. an d cyan with ±la uncerta i nties shown 
as the dashed lin es for| Pcnton, Stocke, fc Shulj (|2004l '). The 
slope obtained bv lPenton, Stocke, fc ShuU (2004 ) is — 1.65=fc 
0.07, w hich is statistical ly different from th at of lLehner et al.l 
(|2007ll (-1.84±0.06) or lDave et all (|200ll ) (-2.04±0.23). As 
we will show in i ]4.4l at least part of the difference likely owes 
to spectral resolution, since higher resolution data yields 
steeper slopes. For this reason it is not straightforward to 
compare to the COS predictions presented here to these 
data, but broadly the slope and amplitude predicted by 
the current simulations are in agreement given the various 
uncertainties. A careful comparison against upcoming COS 
data should provide more discrimination between models. 

At higher column densities, IPenton. Stocke. fc Shulj 

(1200J) found that their data showed a characteristic "dip" 
in the CDD at A'^ei ~ lO^'^cm"^, whic h they point out per - 
sists from high-2 (see their Figure 9). iLehner et al.l (|2007l ) 
likewise saw that the CDD slope is shallower at A^hi > 
lO^^'^cm"^ than for smaller systems. We see a hint of the on- 
set of such a dip in all our simulations, but we cannot accu- 
rately trace it out to Ahi > lO'^^cm"^, since we quickly run 
out of absorbers owing to the rarity of high-A^ni absorbers 
along our randomly-chosen LOS. The origin of this dip is un- 
clear. It occurs at a column density close to that expected 
for overdensities near the boundaries of galaxy halos, which 
leads us to speculate that perhaps these absorbers probe the 
outer parts of galaxy halos where (at least in larger halos) 
some of the gas might be shock-heated to a temperature 
where it would not absorb in Hi. At even higher columns, 
one would then probe through denser halo regions where 
Hi condenses out again, owing to self-shielding, causing a 
"recovery" back to the original CDD. We leave a detailed 
examination of this conjecture for future work. 

The redshift evolution of the CDD shows a steady 
march upwards in line counts at a given A''hi (> lO^^'^cm"^). 
The overall steady evolution reflects the evolution of the 
mean flux decrement over this redshift range, and hence 
is dependent on our procedure of adjusting the ionizing 
background strength to match D{z). The incompleteness 
at A'hi ^ lO^^'^cm"^ becomes more severe at higher-z, as 
line blending becomes more common. Recall that we have 
fixed our resolution and line spread function to the COS far- 
UV channel at all z, which underestimates the data quality 



that is obtainable at 2: ~ 1 (where the COS NUV channel 
has a narrower LSF, albeit with lower sensitivity) and at 
2: ~ 2 (where one can trace the Lya forest in the optical). 
This makes it difficult to compare slopes across a fixed A'hi 
range, but if we restrict ourselves to 10^'^'* — 10^*"cm~^ 
where the CDD is reasonably complete at all z, we ob- 
tain slopes for the vzw model of —1.89, —1.71, and —1.54 at 
z = 0, 1,2, respectively. The CDD is clearly getting steeper 
with time, broadly matching t he observed CDD slope evolu- 
tion f rom high-z (~ —1.5, e .g. iKim. Cristiani. fc D'Odoricd 
l200ll : |janknecht et al.lboOfil 'l to \ow-z. 

The CDDs are remarkably insensitive to the outflow 
model at lower column densities. Wind dynamics, despite 
enriching the IGM and depositing energy over large scales 
in quite different ways, h ave little impact on the diffuse H I 
distributi on. As s hown by Oppenheimer fc Davg (|2008l ) and 
KoUmcicr et al.l ([2003, 120061 ). the typical extent of galac- 
tic winds (~ 100 kpc physical) generally does not extend 
into the diffuse Lya forest. This was noted in the earli- 
est si mulation studi e s of L ya absorption that explored out- 
flows iTheuns et al.l (12003). However, there are signiflcant 
differences for stronger absorbers, as outflows (particularly 
in the vzw case) deposit more cool gas in the outskirts 
of halos, which increases the Lya absorbing cross-section. 
This yields significantly more sub-Lyman Limit systems and 
slightly shallower CDD slopes, and it is also reflected in the 
number counts of strong lines as we will show in H4.3I 



4.2 &-parameter Distributions 

The right panels of Figure [S] show the b-parameter distri- 
butions for our three wind models a.t z = Q, 1,2. The his- 
tograms show the usual Gaussian distribution with an ex- 
tended tail to higher line widths. The median 6-parameter 
for each model is indicated by the tickmark above the curves. 
The typical value is around 40 kms~^ at all z, with a slight 
trend to become smaller at high-z. 

The line widths are si gniflcantly wider t han what was 
found using STIS dat a (e.g. [D ave et al.| l200ll : IWimger et al.l 
(2010 ). Simulations bv lDave et al. (20011 ) predicted a median 



h that depended on A'hi with a typical value of 25 km s~ , in 
good agreement with their STIS data. iLehner et al.l ([20071) 
deter mined a median h = 30 kms~^, also from STIS 
data. 'Paschos et al.l ((2009') used simulated spectra with ap- 
proximately 5 kms~^ resolution and obtained a similar 6- 
parameter distribution. We will show in H4.4l that the larger 
6 values predicted in this work are mostly a result of the 
poorer spectral resolution of COS relative to STIS. Our sim- 
ulated line widths are more compar able to those obtained us- 
ing G HRS (~ 18 kms~^ resolution: IPenton. Stocke. fc Shulj 
[2OOJ). Owing to this sensitivity to instrumental characteris- 
tics, we do not conduct any detailed comparison with data 
here. Nevertheless, we will show in ij5]that line widths still 
contain some information about the underlying temperature 
of the absorbing gas. It may be possible to extract better 
constraints on linewidths using higher-order Lyman lines 
through a curve-of-growth analysis or through simultane - 
ous fitting of multiple transitions fe.g. ILehner et al.l 120071 ') . 
but we leave such a study for the future. 

The line widths, like the CDDs, show almost no sen- 
sitivity to outflow model. This is perhaps more surprising, 
given that the various wind models heat the IGM to dif- 
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Figure 6. Line-of-sight number density of Lya absorbers 
(dN/dz), as a function of redsliift. Top panel shows tlie evolution 
of absorbers with A^hi > lO^^cm"^, while bottom panel shows 
]^q 13.i ^ jY^j ^ 10^^cm ~^. Data from various sources compiled 
bv lWilliger et al.l l|201(Jl are indicated. 



ferent levels (see Figure |4]) , with cw clearly providing sub- 
stantial heating. In ^5.41 we will show that the b parameters 
are only loosely correlated with temperature, reflecting the 
dominance of other sources of line broadening such as Hub- 
ble flow and instrumental resolution. In detail, the cw model 
does have marginally higher median b values of 43.0, 39.4, 
and 38.2 kms~^ at z = 0, 1,2 respectively, as compared to 
vzw which has 42.3, 38.2, and 36.2 kms~^ (almost identical 
to nw), so the expected effect is present but small. We will 
also show in i]4.4l that b parameters are highly sensitive to 
spectral resolution, and hence comparisons of data to mod- 
els (or to other data sets) must carefully account for such 
effects. 



4.3 Absorber Evolution 

The evolution of Lya absorbers provides insights into the 
interplay between cosmic expansion, photoionization, and 
the growth of structure. One simple statistical characteriza- 
tion of this evolution is the redshift-space number density of 
absorbers along the line of sight, dN/dz. 

dN/dz is usually measured wit hin a specified column 
density range. Early data from FOS (jWevmann et al.lll998l ) 
showed a dramatic change in the evolution rate of strong 
(A'hi ^ lO^^'cm"^) absorbers between z > 2 and z <1. How- 
ever, subsequent higher-quality data has muddled the situ- 
ation. Owing to the paucity of such strong absorbers, there 
appears to be substantial cosmic variance for dN/dz among 
the various sightlines. Moreover, more recent data with sub- 
stantially higher resolution generally shows a significantly 
lower dN/dz compared to the FOS results, presumably be- 
cause higher resolution enables more accurate deblending of 
strong lines. It appears that large redshift path length and 
high spectral resolution are both necessary to obtain an ac- 



curate estimate of dN/dz for strong absorbers; COS should 
provide this combination. 

Figure [S] shows the evolution of dN/dz for our three 
outflow models. The top panel shows dN/dz for A^hi > 
10^*cm"^ lines, while the bottom shows 10^^'^ < iVni ^ 
10^*cm~'^ absorbers. Overlaid are a compilation of data 
taken from lWilliger et al.l (|2010n . For strong lines, the points 
across the top of the pane l are the Key Project FOS re- 
sults (|Wevmann et al.lll998l ). The others come from higher- 
resolution data sets, and are all significantly lower than 
the FOS results. In general, the simulations show a rela- 
tively slow decline in dN/dz, with little evidence for any 
kind of "break" in evolution as was se en i n the FOS 
data. Theuns. Leonard, fc EfstathiotJ (| 19981 ) and lDave et al.l 



(Il999f) emphasized that the decline of the U V background at 
z < 2 slows the evolution of the forest relative to a constant 
UVB case. The steadiness o f evol ution in Figure [6] relative 
to Figure 3 of I Dave et al.l (|l999l ) is partly a consequence 
of the narrower redshift range and partly a consequence of 
adopting a UVB that declines more slowly at low redshift 
because of contributions by galaxies. 

The strong absorbers show a clear discrimination be- 
tween outflow models by 2 = 0, with twice as many A^hi > 
10^*cm~'^ absorbers in the vzw model as in the no-wind 
case. The cw model is intermediate. This shows that our 
favoured outflow model is ejecting more gas into the vicin- 
ity of galaxies where strong absorbers arise, and that this 
gas is remaining relatively cool. The difference arises be- 
cause the momentum driven wind model ejects more mate- 
rial from lower-mass galaxies at lower wind speeds, thereby 
providing a reservoir of cooler gas that can cause strong Lya 
absorption. This is qualitatively the same effect that causes 
this model to provide a larger amount of DLA absorption 
and broader DLA kinematics, as described by iHong et al.l 
(2010). The vzw case shows the best agreement with recent 
z ~ data, which is directly tied to its somewhat shallower 
ODD slope, although the uncertainties are still large enough 
to preclude any firm discrimation. COS should yield a good 
observational determination, providing an important diag- 
nostic of galactic outflows. 

The weaker absorbers show much less sensitivity to 
outflows, indicating that winds do not substantially dis- 
turb gas at the densities giving rise to A'hi ^ lO'^^cm"^ 
absorbers. A reasonable flt for all the models is given by 
dN/dz oc (1 + z)'''^, which is significantly s hallower than 
the evolution at high-z fe.g. iDave et al.lll999l '). However, all 
the models generically fail to reproduce the a pparent break 
in evolution at z ~ 1 IjJanknecht et al.ll2006l ). This can be 
traced directly to the evolution of the mean flux decrement 
D (Figure [1}, which the models overpredict around z ~ 1 by 
the same factor of two that they ovthe erpredict dN/dz. In 
fact, the two observational data sets are not unrelated, as D 
is most sensitive to absorbers just below the logarithmic por- 
tion of the curve of growth, namely A'hi ^ 10^^cm~^. If we 
adjust ed the UVB intensity to reproduce the lKirkman et al.l 
(120071 ) D measurements at z = 1 — 1.5, then we would pre- 
dict a break in dN/dz evolution at log(l + z) ~ 0.35. Im- 
proved statistics will thus provide valuable constraints on 
the evolution of the UVB, perhaps demonstra ting significant 
departures from the lHaardt fc Madaul (|200ll ) predictions. 

In summary, the evolution of Lya absorbers will in prin- 
ciple provide strong constraints on the evolution of the ion- 
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izing background (for weaker lines), and perhaps on the na- 
ture of outflows (for stronger lines). However, accurate con- 
straints require larger data sets with good resolution, which 
COS should provide. 

4.4 The Effects of Instrumental Resolution and 

Noise 

Intrumental resolution and noise characteristics of the data 
can have a non-trivial impact on the derived statistics of 
Lya absorbers. This can lead to some confusion when cross- 
comparing different samples, or when comparing simulations 
to data without carefully matching such characteristics. In 
this section we use our simulations to briefly investigate how 
COS-resolution data would compare to an equivalent sam- 
ple of higher STIS-resolution data, and how the Lya forest 
would appear if it were possible to greatly boost the S/N 
of COS data. These are idealized experiments, not readily 
achievable with current instruments, but they illustrate the 
trends one should keep in mind when intercomparing sam- 
ples. 

Figure [7] shows the CDD and &-parameter distributions 
for the vzw model at 2: ~ for our simulated spectra con- 
volved with a COS LSF at S/N=30 and 100, and with 
a 7 kms~^ Gaussian (with 3 kms~^ pixels) intended to 
roughly mimic STIS resolution. The higher S/N or resolution 
both extend the CDDs to lower column densities with an un- 
broken power law, showing that our COS-resolution sample 
with S/N=30 is essentially complete to Nyii = lO^^cm"^. 
Coincidentally, increasing the S/N to 100 at COS resolu- 
tion is roughly equivalent, in terms of CDD completeness, to 
having STIS resolution with S/N=30. At the high-A^'ni end, 
STIS resolution produces quite a few more strong lines. We 
caution that different Voigt profile fitting algorithms can sig- 
nificantly affect the parameters derived for saturated lines, 
so one must be careful in interpreting this plot. However, it 
does indicate that some of the variations seen in dN/dz for 
strong lines may owe to different instrumental characteristics 
for the various data sets. Examination of the true column 
densities in the simulated absorbers shows that the high col- 
umn density systems identified in the STIS-resolution spec- 
tra are generally real; at COS resolution AuTOVP mischar- 
acterizes them as lower column density systems with larger 
6-parameters (and thus less saturation). Finally, the slope 
of the CDD at lO" < JVhi < 10^^-^cm"^ is steeper in the 
case of higher S/N or resolution; for the former, we obtain 
a slope of —1.73 (as opposed to —1.70 for the original vzw 
case), and for the latter we get —1.82. Hence, high resolu- 
tion produces a significantly steeper CDD slope, wh i ch can 
help explain why the CDD slopes in iDave fc Tripd (|200ll ') 
are steeper. 

The 6-parameter distribution shows strong sensitivity 
to spectral resolution, but increasing S/N does not have a 
major impact on b parameters. At STIS resolution, the me- 
dian &-parameter is below 30 kms~ ^, and the d istribu tion 
roughly agrees with that found in iDave et al.l (|2001h . At 
COS resolution, lines are significantly broader, resulting in 
a median greater than 40 kms^^ and a much more promi- 
nent tail to high 6 values. The cyan dot-dashed curve in 
Figure [7] shows the result of taking each line in the STIS- 
resolution distribution and increasing its 6-parameter by 
adding ^2(17^ — 7^) ~ 22 kms~^ in quadrature to account 
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Figure 7. Effects of signal-to-noise ratio and spectral resolution 
on Voigt profile parameter distributions in our fiducial vzw run at 
z = 0. COS resolution assumes the COS LSF in the FUV channel, 
which is Si 17 kms~^ Gaussian width (lo"), while STIS resolution 
assumes a Gaussian with 7 kms~^ width. Top panel shows col- 
umn density distributions multiplied by A'^pjl i demonstrating that 
the turnover at low-Afjjl a-t lO^^cm"^ owes to a lack of sensitivity; 
either increasing S/N or resolution allows one to probe further. 
There are also differences at the high-Afjji end. Bottom panel 
shows the 6-parameter distributions; these are not very sensitive 
to S/N variations, but are dramatically sensitive to spectral reso- 
lution. Here we also show the result of convolving the linewidths 
identified in the STIS-resolution spectra with a 22 km/s Gaussian 
to illustrate the effect of moving from STIS to COS resolution; 
this demonstrates that much of the line broadening in the COS- 
resolution case owes to its coarser spectral resolution. 



for the difference in resolution between COS and STIS. The 
peak of the distribution moves much closer to that of the 
COS distributions, suggesting that for many lines the larger 
COS 6-parameter is a simple consequence of convolving the 
line with a broader LSF. However, the COS distributions 
still have many more lines at 6 J5 60 kms~^, showing that 
these broad lines must arise mainly from blends that cannot 
be resolved by the COS FUV channel. Because of the effects 
of line-blending, we think it is generally better to test mod- 
els by applying the LSF to simulated spectra rather than 
attempting to correct measured 6-parameters on a line-by- 
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line basis. The exception might be in cases where higher or- 
der Lyman series lines provide additional constraints on the 
fe-parameters via curve-of-growth analysis. More generally, 
this illustrates that the b parameters are particularly sen- 
sitive to spectral resolution, and any comparisons between 
data sets, or between data and models, should take special 
care to match the instrumental line spread function. 



5 Lya ABSORBER PHYSICAL CONDITIONS 

We now examine the relationship between absorber proper- 
ties and the physical conditions of the absorbing gas, namely 
its temperature, density, and ionization state. The rela- 
tively simple nature of Lya forest absorption, particularly 
for weaker absorbers, results in fairly tight relationships that 
offer insights into the evolution of baryons in the IGM. 

To assign densities and temperatures to individual ab- 
sorbers, we find the location of peak absorption and assign 
to the absorber the H l-weighted p and T of the nearest 
pixel. This is not a perfect procedure, as, e.g., the width 
of the line might be determined by coincident hotter gas 
superposed on stronger absorption from colder gas. But it 
gives a reasonable idea of the physical conditions present 
in the absorbing gas. To better isolate real trends, we only 
consider absorbers whose formal x^ errors for A'^ei and b as 
reported by AuTOVP are fractionally less than 20%, which 
mostly removes very weak absorbers whose parameters are 
not well-determined. For studying linewidths, we include 
only absorbers with A'^ei > 10^'' cm^^, for which we have 
a relatively complete sample. We focus on the momentum- 
driven wind case and point out where substantial differences 
arise in the other wind models. 



5.1 Absorbers in Phase Space 

Figure [8] shows how Lya absorbers populate and trace cos- 
mic phase space (here shown as temperature vs. hydrogen 
number density). In each panel except the upper left, our 
simulated absorbers between z = — 0.2 are indicated as 
circles, placed using their Hi-weighted densities and tem- 
peratures. Solid lines demarcate the phase boundaries as 
defined in i]3.2l and shown in Figure (3] and the vertical dot- 
ted line indicates the cosmic mean density of hydrogen. Our 
absorber sample overlaps the shading, which represents the 
simulated distribution of Hi. However, because random sight 
lines sample the Lya forest in a volume-weighted fashion, 
there are more weak Lya absorbers tracing lower densities 
even though most of the cosmic H I resides in rare, strong ab- 
sorbers tracing the condensed phase. Conversely, absorbers 
are not detectable at densities below uh = lO^'' cm~'^, 
where the vast majority of the IGM volume exists, and this 
detectability threshold is the reason that the low-z Lya for- 
est is observed to be sparsely populated. 

The upper left panel displays the cosmic distribution of 
baryons in phase space, color-coded by ionization fraction. 
The remaining 5 plots contain up to three additional dimen- 
sions of information within two-dimensional cosmic phase 
space: Shading for global Hi content (i.e. the product of 
the intensity and color maps in the top left panel), circle 
size scaled to absorber A^hi, and circle color for the various 
quantities, as follows: Neutral fraction (upper right), total 



hydrogen column (middle left), absorber path length (mid- 
dle right), A'hi (lower left), and linewidth (lower right). It is 
clear that there are strong correlations amongst the various 
physical and observational quantities. In subsequent sections 
we will quantify some of these trends in more detail and at- 
tempt to connect observables to physical quantities within 
this plane. 

The Hi neutral fraction in the upper right panel is 
computed for each absorber as an optical depth-weighted 
neutral fraction from all SPH particles that contribute to 
the absorption. Note that we assume ionization equilib- 
rium; this is expected to be a very good approximation 
in the Lya forest. It is clear that, as expected, the IGM 
is highly ionized, with the weakest absorbers showing neu- 
tral fractions of ~ 10~® cm~^. Neutral fractions increase 
towards lower temperature and higher density, such that 
absorbers at ~ lO^^cm"^ (the typical absorber at the cross- 
ing point of the phase division lines) show a neutral frac- 
tion of ~ 10~^ cm~^. The neutral fraction increases rapidly 
along the condensed phase, leading quickly to very strong 
absorbers. 

The middle left panel of Figure [5] shows the total hy- 
drogen column density Nh, computed by dividing A'hi by 
the neutral fraction. A remarkable trend emerges that the 
vast majority of Lya forest absorbers have Nh between 
10^*cm~^ and lO^'^cm"^, except for the hottest or densest 
absorbers. This trend was noted observationally in a small 
sample bv lProchaska et al.l (|2004l ). which they pointed out 
was in general ag reement with simulation predictions from 
I Dave et al.l (|200l|). Here the trends of Nh versus overdensity 
can be more clearly seen; it is double-valued in the sense 
that Nh is lower at both low and high overdensities, and 
it peaks at intermediate overdensities, particularly at high 
temperatures. 

The middle right panel shows an estimated size of Lya 
absorbers, obtained by dividing Nh by the physical density. 
Weak absorbers (~ lO^^cm"^) have sizes up to ~ 1 Mpc, 
and this quickly drops to ~ 100 kpc for ~ 10^*cm~^. This 
is broadly consistent with the large coherence l engths ob- 
served in quasar pair sightlines (e.g. iDinshaw e t al. 1993; 
ICasev et all 120081 ') , though the elongated geometry of ab- 
sorbing structures also affects apparent coherence l engths. 
The s izes are consistent with analytic estimates bv ISchavd 
(J20i)l') from assuming that Jeans smoothing establish es the 
structure of the low-density IGM. (See Peeples et al.i [20101 ] 
for evaluation of this hypothesis against simulation results, 
albeit at higher redshift.) The absorber lengths drop rapidly 
with increasing overdensity, such that strongest absorbers in 
the condensed phase have sizes ~ 1 kpc or less. This explains 
why absorbers from this phase are rarely seen despite the 
substantial reservoir of neutral gas in the condensed phase 
— their cross-section is remarkably small. 

The two lower panels consider absorber observables A'hi 
and b. The lowest column density systems lie along the low- 
density, high-temperature envelope of detectable Hi. The 
trend of increasing A'hi towards low temperatures and high 
densities mimics that in the ionization fraction, such that the 
typical absorber at the bound/unbound division has A^hi ~ 
lO^'^cm"^ (and lower at higher temperatures). The strongest 
absorbers generally arise in condensed phase gas, though 
they are rare in these random LOS. 

The 6-parameters in the lower right panel also show 
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Figure 8. Cosmic pliase diagrams for tlie vzw run a.t z = 0.1. Upper left sliows all baryons, colored by Hi fraction. Other panels are 
shaded by flHly binned in 0.1 X 0.1 dex pixels. Solid lines demarcate the four cosmic phases as in Figure|31 while the dotted line indicates 
the cosmic mean density. In each panel, the circles show the locations of z = — 0.2 absorbers in phase space, with the area of the circle 
proportional to A^hi (see lower left panel for scale). The color scale in each panel corresponds to different quantities: Ionization fraction 
(upper right), total hydrogen column density Nh (middle left), absorber size (middle right), A^hi (lower left) and linewidth (lower right). 
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Figure 9. Baryonic physical density in units of the cosmic 
mean (p/p) versus column density Afjji for absorbers from the 
momentum-driven wind model, at z = (top), 2 = 1 (middle), 
and 2 = 2 (bottom). In each panel, the solid line shows the re- 
lation from equation |3] at 2 = (black), 2 = 1 (magenta) and 
2 = 2 (cyan). Absorbers are color-coded as diffuse (blue), WHIM 
(green), hot (red), or condensed (cyan). Point sizes are scaled oc b. 



trends in phase space, but they are not as distinct as with 
column densities. Narrower lines tend to arise in lower- 
density and lower-temperature gas, but there is a significant 
scatter. Wide lines can occur in hot gas, but such gas can 
also contain narrower lines. Wide lines can also occur in very 
strong absorbers that are heavily saturated, arising in con- 
densed gas. Overall, b parameters trace temperatures fairly 
loosely, a point we will quantify further in i ]5.4l 

Overall, Lya absorption along random lines of sight can 
trace present-day baryons from close to the mean overden- 
sity up to overdensities of ~ 10* or more, and to tempera- 
tures well above 10^ K. There are strong trends of ionization 
fraction and absorber size, which together with the underly- 
ing baryon distribution in phase space yield a strong trend 
of physical parameters versus column density, and a weaker 
but still noticeable trend versus linewidths. In the next two 
sections we explore these trends more quantitatively, to bet- 
ter understand how observables obtained from COS spectra 
can reveal the physical conditions of the absorbing gas. 



5.2 Column density versus physical density 

The evolution in the relationship between column density 
A^Hi and overdensity 5 is governed by the interplay between 
cosmic expansion, which lowers the physical density (at a 
given S) and thus lowers the recombination rate, and evolu- 
tion of the UV background, which drops at lower redshifts 
and thus lowers the photoionization rate. I Dave et al.l (|l999l ) 
showed that much of the observed evolution in the Lya for- 
est can be understood through the evolution in this relation, 
because this interplay causes a given strength absorber to 
correspond to a higher overdensity at lower redshift. Compli- 
cations arise because, by z; = 0, a significant amount of IGM 
gas is heated coUisionally by shocks rather than b y photoion- 
ization, resulting in departures from the FGPA ( Zhan et al.l 
2005 ). Nevertheless, this A^hi — S relation is central to un- 
derstanding the physics of the Lya forest in a structure for- 
mation context. 

Figure |9] shows that a fairly tight A''hi — 5 relation is es- 
tablished at all redshifts, as found by previous studies. The 
three panels show this relation for the momentum-driven 
wind model at z = 0, 1, 2. Absorbers are categorized by cos- 
mic phase, with a color code following Figure [J] namely dif- 
fuse (blue), WHIM (green), hot halo (red), and condensed 
(cyan). The lower envelope reflects gas that is purely pho- 
toionized, as shock heating will move absorbers to lower A^hi 
for a given overdensity by lowering the neutral fraction. 

The diffuse-phase absorbers follow a tight relation in 
overdensity versus temperatur e owing to photoionization 
heati ng and cosmic expansion (JHui fc Gnedinll 19971 : ISchavd 
[2OOIJ). We perform a least-squares fit to the diffuse absorbers 
as a function of column density and redshift (from z — 0—2), 
restricting ourselves to absorbers with T < 10*'^ K to re- 
move mildly shocked systems, and obtain the following re- 
lation: 



p 35.5 ±0.3 
P 



/?■ 



V10i4cm-2/ 



0.741±0.003 



10 



(-0.365±0.009)z 



(3) 

Recall that fr is the factor by which optical depths 
have been multiplied r elative to those computed for the 
iHaardt fc Madaul (|200ll ) ionizing background, to match the 
mean flux decrement; it is 2/3 for our momentum-driven 
wind model. The quoted uncertainties are on means, and 
they do not reflect the spread amongst the individual ab- 
sorbers. This fit is shown at each redshift in each panel, to 
help visualize the rate of evolution. Our slop e for p vs. Ah i 
is similar to the slope of two-thirds derived bv lSchavd (|200ll ) 
from Jeans smoothing arguments. Furthermore, the redshift 
evolution is in excellent agreement with predictions from his 
analytic model. 

The WHIM absorbers generally have somewhat higher 
density at a given A^hi, since the gas has been shock-heated; 
this trend is exacerbated for the hot halo absorbers. As a 
result, at overdensities above 10 or so, corresponding to 
N-Ri <; lO^^cm"^ at z ~ 0, there begin to be significant 
departures from equation ^. The departures from this re- 
lation also become more prominent with time, since more 
IGM gas becomes shock-heated. Condensed phase gas typi- 
cally gives rise to the strongest Lya absorption systems seen 
at the highest overdensities. At z = 2 there is a shelf at the 
lowest densities, largely because the underdense absorbers 
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Figure 10. Summed column densities per unit path length as a 
function of column density, subdivided by absorber phase. The 
top panel shows absorbers in the range z = — 0.5, and the 
bottom panel shows z = 1.5 — 2. 



have larger velocity widths and are thus more difficult to 
identify at fixed A^hi and S/N. 

Equ ation ([3]) is s i milar in fo rm to the (by-eye) fits pre- 
sented in lDave et al] (|l999l ') and iDave eT^ (|200ll '). Differ- 
ences arise because the ionizing background strength em- 
ployed here is slightly different, and here we consider only 
the diffuse absorbers. We have checked that the above equa- 
tion holds reasonably well for the diffuse absorbers in the 
no wind and constant wind cases (when the /t factor is in- 
cluded), since these absorbers are essentially unaffected by 
outfiows. The other wind models are qualitatively similar in 
terms of their other phases as well, though in detail there 
are minor differences that we will discuss in the following 
sections. The insensitivity to winds further emphasizes the 
fundamental nature of the p — Nui relationship in describing 
the physics of the Lya forest. 



5.3 Absorption phase 

We now examine the cosmic phase of Lya absorption as a 
function of column density. Figure \TU\ shows a sum over H I 
column density for lines in bins of A'^hi for our momentum- 
driven wind model. At both low {z — 0.25) and high 
{z = 1.75) redshifts, there is a clear transition in which 
total Hi absorption is dominated by diffuse absorbers at 
low A^Hi and by condensed absorbers at high Ahi- The 
transition occurs at A'^ei ~ lO^^cm"^ for z = 0.25 and 
^ 10^''-^cm-2 for z = 1.75. WHIM and hot halo absorbers 
are sub-dominant at all column densities, though they are 
much more common at z = 0.25 than at z — 1.75. While the 
column density demarcation is less clear for these phases 
than for diffuse/condensed, WHIM absorbers typically have 
A'hi < lO^^cm"^, while hot halo absorbers sometimes have 
larger A'^hi- 

While the broad trends are similar for other wind mod- 



els, in detail there are substantial differences particularly for 
absorption in bound phases. If we consider absorbers from 
10^* < Nm < lO^'^-^cm"^ at z = - 0.5, then the cw model 
contains only half the absorption, and nw only ~ 40% of 
the absorption, relative to vzw. This reiterates the result in 
Figure[B]that high-A'Hi absorbers provide the most sensitive 
probe of galactic outfiows in the Lya forest. 

In summary, the low-z; Lya forest contains absorbers 
from all cosmic phases, but diffuse absorbers dominate at 
A'hi ^ lO^^cm"^ and condensed absorbers dominate at 
higher A'^ei- WHIM and hot phase absorbers are present, 
but highly subdominant in total absorption. Next, we assess 
whether linewidths add enough information to pinpoint the 
absorbers arising in WHIM and hot halo gas. 

5.4 Linewidths and temperatures 

In the absence of other sources of line broadening, the b- 
parameters will refiect the temperature of the underlying 
absorbing gas. A particularly promising application of this is 
tracing the missing baryons with wide Lya lines because, un- 
like relying on high- ionization metal lines su ch as Ovi (e.g. 



Tripp. Savage, fc JenkinjIJOOol ) or O vii (e.g. iNicastro et al.1 
[5|), Hi linewidths are independent of metallicity. How- 



200, 



ever, in prac tice, linewidths are no t dominated by thermal 
broadening ([Weinberg et al.l 1 1997c ), although the t hermal 
contribution is higher at low-2 ( Dave fc Tripdl200ll ). Here 
we investigate the relation between linewidths and temper- 
ature in our simulations, to better understand how robustly 
one can trace the WHIM, and hotter IGM gas in general, 
using wide Lya lines. 

Figure 1111 top panel, shows the distribution of 
absorbers' 6-parameters versus temperature from our 
momentum-driven wind simulation at z ~ (other wind 
models are similar). Absorbers are color-coded by cosmic 
phase as in Figure U Recall that we do not fit a continuum 
to the simulated spectra, hence there is no additional uncer- 
tainty in this regard in the identification of broad absorbers 
in the models. However, such uncertainties may be impor- 
tant when analysing observations, as broad absorbers may 
be partially "fit out" when fitting a continuum to unnormal- 
ized data. There is a definite trend for wider absorbers to 
arise in hotter gas, as shown by the running median (ma- 
genta curve), but there is also substantial scatter as indi- 
cated by error bars which show la dispersion values about 
the median. Hence, a wide Lya line is more likely to arise 
in hotter gas, but it does not always do so. 

If lines were purely thermally broadened, then the H I 
lincwidth would be related to the temperature by b = 
12.8^/T/10''K kms"^. This relation is shown as the solid 
diagonal line. A small fraction of WHIM and hot absorbers 
are narrower than allowed by thermal broadening, either 
because of 6-parameter uncertainties in the Voigt profile 
fits or because of uncertainties in assigning temperatures 
to absorbers as described above. Generally, however, most 
absorbers' widths significantly exceed the thermal broaden- 
ing value. The running median shows a similar trend as the 
thermal broadening track, but it is higher by a factor of 
1.5 — 3 in 6, until a turnover at the very largest b values. 
For t he narrow diffu s e abso rbers, we obtain the same result 
as in iDave fc Tripd (|200ll ) that thermal broadening con- 
tributes equally to other sources (in quadrature) in setting 
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Figure 11. (Top) Temperature vs. fe-parameter for Lya ab- 
sorbers, color-coded by absorber phase, from the vzw simulation 
a.t z = 0. Symbol size is scaled by ^/A^jji/lO^*. Diagonal solid 
line shows the relation for pure thermal broadening; some ab- 
sorbers are narrower than "allowed" owing to uncertainties in 
Voigt profile fitting and the assigning of temperatures to individ- 
ual absorbers. Histograms in b and T are shown along the bottom 
and side, respectively, divided by absorber phase, with tick marks 
above the curves identifying the median values. For visibility of 
the other phases, the diffuse-phase histograms have been reduced 
by 5x; otherwise, it would dominate at essentially all b values. 
Magenta curve shows a running median, with Ict dispersion, for 
the entire sample. (Bottom) Like top, but showing Hi column 
density in place of temperature. 



the linewidths. For wider absorbers, the thermal contribu- 
tion becomes smaller. In general, the dominant contribution 
to the 6-parameters is not thermal broadening, but rather 
Hubble broadening in the expanding gas or line blending 
owing to limited spectral resolution as discussed in H4.4I 

Along the bottom and side of Figure [TT] are displayed 
absorber histograms in b and T for each cosmic phase, with 
the median values indicated by tick marks. The diffuse-phase 
histogram has been multiplied by 1/5 along each axis for 
visibility. The ones along the temperature axis show the ex- 
pected behavior, with hotter phases having higher median 



temperatures. The increasing T with b is broadly reflected in 
the histograms along the &-axis as well, with WHIM and hot 
absobers being at higher b values than diffuse and condensed 
ones. However, the scatter between b and T is quite signifi- 
cant; the histograms of the various phases span a wide range 
in b and have a substantial overlap. Moreover, the diffuse ab- 
sorbers essentially dominate at all b values (remember that 
this histogram has lowered by 5x) simply because they are 
so numerous, despite the fact that their median linewidth is 
substantially smaller than all the other phases. Hence, it is 
unfortunately not straightforward to assign a gas tempera- 
ture purely from the b parameter. A key implication is that 
searches for WHIM gas with COS by selecting only broad 
Lya lines is likely to pick up a significant fraction of diffuse 
absorbers, and care must be taken to avoid double-counting 
baryons in such studies. 

Both the column density and linewidth provide some 
guidance in determining the cosmic phase of the absorbing 
gas. It stands to reason then that combining the two observ- 
ables may provide more accuracy in identifying the cosmic 
phase. The bottom panel of Figure [TT] shows the absorbers 
in A^Hi — b, color-coded by absorber phase. Broad, weak ab- 
sorbers with 6 > 60 kms~^ and A^'hi < lO^^cm"^ include 
a relatively large fraction of WHIM absorbers. However, a 
significant fraction of absorbers in this range are still dif- 
fuse, and many WHIM or hot halo absorbers have lower 
6-parameters than this threshold. 

In summary, the linewidths trace temperature only 
fairly loosely, and for any individual absorber it is difficult 
to pinpoint the cosmic phase purely from observables. Nev- 
ertheless, one can assign phases probabilistically based on 
their location (NHi,b), at the penalty of introducing some 
model dependence. With the large ensembles of absorbers 
expected to be obtained with COS, this approach can pro- 
vide some constraints on the content of all baryonic phases 
as traced by Lya absorption. 



6 CONCLUSIONS 

We study the low-redshift IGM as traced by Lya absorbers 
in 48/i~'^Mpc, 2 x 384'^-particle cosmological hydrodynamic 
simulations with three prescriptions for galactic outflows. 
This w ork builds on ou r previou s works o f a de cade ago, 
namely iDave et all lfl999 'l and Dave et al.1 (1200 ll ). by em- 
ploying simulations with 216 x more particles and 80 x more 
volume, using improved numerical methodology, and most 
significantly, testing the impact of large-scale galactic out- 
flows on the IGM as traced by Lya absorption. Our scenario 
for the nature and evolution of Lya forest absorbers from 
z — 2 ^ remains basically unchanged: most absorbers 
arise in low-overdensity, highly p hotoionized gas that traces 
uncoUapsed large-scale filaments. iDave et al . (1999) argued 
that much of the observed evolution in the Lya forest can be 
understood in terms of the evolving relationship between H I 
column density and the overdensity of the absorbing gas; our 
latest simulations yield a similar slope and evolution of this 
relation. Outflows are found to have a small impact on the 
diffuse IGM, but they do have a signiflcant effect on bound- 
phase absorbers that arise nearer to galaxies. Our statistical 
results for Lya forest absorbers are obt ained by applying 
the AutoVP line-fitter (JDave et al.lll997^ to artificial spec- 
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tra with a resolution and S/N similar to those expected for 
typical COS observations. 

Our main results are as follows: 

• Baryons at 2; ~ are divided comparably amongst 
the diffuse phase, the WHIM, and the bound phases (con- 
densed, hot halo, and stars). The WHIM here is defined as 
low-density (outside of halos) gas with T > 10 K, a no- 
table change from previous definitions that impose a maxi- 
mum temperature but not a maximum density. In detail, our 
favoured momentum-driven outflow model has 41% in dif- 
fuse gas, 35% in bound phases, and 24% in WHIM at 2 = 0, 
with values for other wind models and redshifts shown in 
Figure 2] 

• Our favoured momentum-driven wind (vzw) model 
shows similar evolution in diffuse (unbound) IGM baryonic 
phases to the no- wind (nw) case. In contrast, the condensed, 
hot halo, and stellar fractions are all noticeably different be- 
tween those two, and between either of these models and 
the constant wind (cw) model. Furthermore, cw significantly 
heats the IGM, resulting in more WHIM and less condensed 
and diffuse gas. 

• The vzw model at COS resolution shows a column den- 
sity distribution with a power-law slope of —1.70. The slope 
is slightly steeper for nw and cw, and it is steeper for higher- 
resolution spectra. While formally all slopes are within ~ 2a 
of each other, recall that the lines of sight are identical 
through each simulation, so the differences are fairly robust. 

• The shallower vzw slope arises because it produces more 
strong (A^Hi > lO^'^cm^^) lines at z = 0, which agrees better 
with observations of dN /dz evolution. Weaker lines show no 
differences in dN/dz among wind models, but they do indi- 
cate a di screpancy around z ~ 1 t hat may reflect deviations 
from the lHaardt k, Madarj ((20011) ionizing background evo- 
lution. 

• The linewidth (6) distribution shows a median around 
~ 40 kms~^ at COS resolution, and it is insensitive to the 
wind model. It is highly sensitive to spectral resolution, even 
at large b values, and much of the broadening owes to in- 
strumental resolution (note that our 6-parameters are not 
corrected for this). 

• Diffuse Lya absorbers follow a relation between TVhi 



and overdensity given by p/ p oc N-^i 10 



in general 



agreement with our previous findings. This results in a re- 
gion of a given overdensity giving rise to a lower column 
density systems at higher redshift. Departures from this rela- 
tion become more prominent at low-z owing to an increased 
amount of shock heating. 

• The bulk of Lya absorbers have total hydrogen columns 



of 



10^ 



cm , refiecting the trend that the Iow-A'hi 



absorbers are more ionized (neutral fraction ~ 10~ cm~ ) 
and larger (hundreds of kpc), while high-A^ni absorbers are 
less ionized (neutral fraction ~ 10~^ cm~^) and smaller (few 
kpc). 

• Linewidths only loosely reflect the temperature of the 
absorbing gas, because typically linewidths are 1.5 — 3 times 
the thermal broadening value. Hence, wide lines do not ex- 
clusively track WHIM or hot gas, though trends in p and T 
with both h and A^hi could be used to assign absorbers to a 
particular cosmic phase in a statistical sense. 

• Improved dN/dz statistics for weak lines (A'^ei < 
10^*cm~^) will improve determinations of the UV back- 



ground intensity at 2: < 2. Improved statistics at higher 
column densities can provide strong tests of galactic outflow 
models. 

The installation of COS on Hubble promises a great leap 
forward in our understanding of low-redshift Lya absorp- 
tion and the IGM from which it arises. Besides testing the 
basic paradigm of the cosmic web origin for weak Lya for- 
est absorbers, COS Lya data will provide constraints on 
the baryon content in various cosmic phases, ionizing back- 
ground evolution, and the nature of galactic outflows. Ad- 
ditional constraints will be provided by the distribution of 
IGM metals that can be probed to much greater precision 
with COS (the subject of Paper II in this series) and by the 
relation between galaxies and absorption systems (covered in 
Paper III). The combination of these constraints with over- 
all Lya absorber statistics will provide powerful diagnostics 
for the mechanisms of galactic outflows. Comparisons of the- 
ory and observations in the coming years will critically test 
our prevailing paradigms and prescriptions, thereby provid- 
ing valuable insights into the evolution of the majority of 
baryons over cosmic time. 
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